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Abstract 

In this paper we consider Cronin momentum behaviour for p + A, d + A and A + A colh- 
sions in saturation model. Our analysis shows that Cronin momentum behavior at different 
rapidities and energies, can be related with scaling law using simple dimensional considera- 
tion. Using exact numerical solution of Balitsky-Kovchegov equation we show that although 
this dependence is slightly different for McLerran-Venugopalan and Balitsky-Kovchegov 
definition of gluon distribution function in simple model in this case dependencies is al- 
most the same (i.e ratio of Cronin momentum calculated using these gluon distribution 
functions is big constant). This can be used to experimentally distinguish this two variant 
of gluon distribution function definition in saturation model and choose the right one. 



1 Introduction 

In our previous paper ^ it was found that in saturation model Cronin momentum (i.e. 
momentum at which Cronin ratio 2 have maximum) has simple behavior in central rapidity 
region which holds in p + A and ^ -|- A collisions 

ln{qc) = a + bln{^) (1) 

where parameters a and b are approximately defined as: 

a=^^^^^HA)+lniql,), (2) 

1 



^ 0.1304. 



2 + A 

with geometric scaling parameter A = 0.3. It was also showed that for different defini- 
tions H9llO() of unintegrated gluon distribution function slope b is slightly different. 

Relation (Q) was inspired by simple expression Q for saturations scale Qs{x) based on 
geometric scaling effect and dimensional considerations Q. 

logg2(y) = XY + log{Q^,iYo)x^) (3) 

Ql = Pqc (4) 

Since our calculation was done using certain model for dipole forward scattering amplitude 
and approximate expression for cross section proposed in \3[ then it is most probably that 
does not holds if dipole forward scattering amplitude defined with solid theoretical 
background will be taken into account. It could be shown below that this indeed true and 
simple formula Q will be asymptotically restored only in high energy region. 

There is mostly the only one way to define dipole forward scattering amplitude: 
Balitsky-Kovchegov evolution equation |4|, which describes the evolution of the dipole 
forward scattering amplitude N{r,y) of a QCD dipole of transverse size |r| with rapidity 
Y = ln{l/x), 

dNi\r\,Y) ^ («,iV,) /■^2„ r2 
dy 27r2 

here for simplicity we suppose that nucleus is cylindrical (i.e. we completely drop impact 
parameter b dependence from forward dipole scattering amplitude) and strong coupling 
constant is fixed. There is also many its derivatives but we'll leave them for now since 
most of them aren't established very well and bring to additional complexity in numerical 
calculation. It was shown in series of papers |6| that equation ^ belongs to Fisher- 
Kolmogorov-Petrovsky-Piscounov class [7j if N{r, x) transformed to momentum space. It 
can be easily shown |S] that instead of simple scaling law (jSJ we will have the following 
one: 



/ d'z^^:^(Ar(|r-z|) + iV(|z|)-iV(|r|)-iV(|r-z|)iV(|z|)) (5) 



,2,..^_„-X(7c),, 3 3 / 2vr 1 



logQ (y) = a^Y - logy - ^^^^ + 0(l/y) (6) 
7c 7^ V ax (7c) VY 

where xil) = 2'i/'(l) — ^"(7) ~ V'(l ~ 7) is well known BFKL kernel, and 7c = 0.6275... is 
the solution of equation x(7c) = 7cX(7c)- It is oblivious that equation does not holds 
in this case but probably will be restored at high energy. 



2 



It is well known that in saturation model there is two ways to define unintegrated 
gluon distribution function based on N{r,Y): Balitsky-Kovchegov way and McLerran- 
Venugopalan way: 



as(27r) 



l') = / e-^'- VlNG{r, x)) (7) 



<^T{-,<l') = / d're-^^^- -,Na{r,x)) (8) 

as{2-Ky J 

or using rotation invariance of N(r,x) we can derive 



4>A{x,q^) = f^^J^2 L drJo{qr)rVtNG{r,x) (9) 







2, ASaCf r w _ 



= — — 2 / drJo{qr)-NG{r,x) (10) 

This two definition gives similar momentum dependency for unintegrated gluon distribution 
function at high momentum (i.e pure bremsstrahlung) but has a quite different behaviour 
at low momentum (fig. 

We could expect that these different definitions for gluon distribution functions with 
N{r, x) defined throw lead to same dependence of Cronin momentum. Since equation 
was written with fixed strong coupling in mind it is not oblivious why we even should have 
different dependency for Cronin momentum from H9I10() . However since initial condition 
already contains soft scale Aqcd to provide right ultraviolet asymptotic for unintegrated 
gluon distribution function and l[9|). (|10jl have very different behaviour at infrared momenta 
(in case of @ unintegrated gluon distribution function approaches to while q ^ 0, but 
with (fTn|) unintegrated gluon distribution function approaches to oo) it is clear that using 
© and (|10j) will lead to different behaviour for unintegrated gluon distribution function. 

In saturation framework gluon production cross-section in AB collisions can be written 
in simple fcf-factorized form: 

"^'"^'^ - ^"^^ [d'kMxi,Q')Mx2,{q-kf) (11) 



cPq dy Cp q^ 

where (j)A.B is unintegrated gluon distribution of nucleus and proton and xi, X2 are defined 
by equation 

xi = ^e-^ X2 = ^ey. (12) 
Of cause in leading log approximation (|lll) reduces to the following form 

da^B 2as 1 ( . (^^^^q2^GB{x2,q^)+GA{xr,q^)<l,B{x2,q^)) (13) 



d^q dy Cp q^ 
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Figure 1: Unintegrated gluon distribution function defined by Balitsky-Kovcfiegov(red) and 
McLerr an- Venugopalan (green) 
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where G{x,q^) is gluon distribution function. But it is clearly seen that ()13|) and have 
very different behaviour at low momentum when unintegrated gluon distribution have the 
form Q, and usage of should be avoided for momentum dependency calculation. 
Having at hand cross section for A + B scattering Cronin ratio can be easy written 



There are three components we have for numerical calculation: cross section (jllj) . uninte- 
grated gluon distribution functions (|9|l. l|10|) . forward scattering amplitude N{r,x) evolved 
by Balitsky-Kovchegov equation ©• We address them in order. 

First lets consider cross section calculation. There is a few problems with direct nu- 
merical calculation of integral in cross-section It should be stressed that even if 
unintegrated gluon distribution function defined in © has no singularity at = 0, ((Tm) 
does. So we should write numerical integration procedure which symmetrically avoids both 
singularities in (|1H) . one from each gluon distribution function. Since 4>{k) oc p- at high 
momentum k it is hard to calculate integral in directly due slow downfall at high 
momenta. Of cause it is not a problem in analytical calculation since this integral is a good 
defined one. Since direct computation has have before-mentioned problems it is clear we 
should use some other method. Moreover even if we decide to do direct computation for 
this integral it will cause unnecessary loss of precision. 

It is easy to calculate integral in by the following method. Lets transform integra- 
tion area from to something more foreseeable, like (numerical integration in infinite 
boundaries is complicated procedure by itself and should be avoided if possible). It will be 
convenient if resulting integrand will have same order of magnitude in all and we chose 
transformation which satisfy this condition. Lets divide singularities i.e. separate integral 
into two pieces. It can easily be done with following prescriptions 



Rab 




(14) 



With this Cronin momentum can be easily calculated. 



2 Numerical Methods 



da^^ _2as 1 f 

(Pq dy Cf Jc^xqap 



d^k(j)A{xi,q^) (I)b{x2, (q - kf) + 



(15) 
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where go is some value with property < qq < p (for simphcity we'll take go = f- Now 
we'll make the following parameter transformation in first term in ()16() 



(16) 



where ipQ = arcos{'min{l, {qo/q)) Similar transformation should be done in second therm. 
It is clearly seen that equation ()16() transform vector q into pair (x, y) which belongs to . 
Of cause this transformation gives additional measure to integrand in (jlbj) which equals to 
(vr — ipo) ■ After this transformation both resulting integrals could be calculated very 
fast by adaptive Cubo method. 

Lets consider unintegrated gluon distribution function now. It clearly seen that z in- 
tegration in Q can be easily done since V'^ efficiently cuts off integration area. However 
it is not the case for (|lflj) since integrand falls slowly with increasing of z. One of possible 
method for numerical calculation of ()1U() is just cut off z integration at some big z and 
analytically calculate the rest. This however leads to quite complicated formula for tail 
term. We choose here more simpler (from numerical point of view) method. Since roots of 
Bessel function Jq{x) is well known and approximately equals to 

x, = {i + ^)T: (17) 

we can define set of finite integrals from to Xj+i 

= -nr^ / drUqr)-NG{r,x) (18) 
as{2-KY Jo r 

It could be easily showed that Si have monotonic behaviour at high i. By applying Eck- 

ler process to pair {j,Si) and using polynomial extrapolation we can calculate (jlUj) with 

any accuracy. It should be noted that using more precise value for Bessel function roots 

gives substantial speedup for this procedure and therefore we'll use more precise McMahon 

formula which in next order gives. 

^^ = ^^ + 8^--3^' ^''^ 
And the last and most important component is forward gluon scattering amplitude 
N{r, x). It should be noted what Balitsky-Kovchegov equation can be solved much simpler 
in momentum space rather coordinate space due its analytical properties. However we'll 
solve it in coordinate space. And we have solid reason for this. Developed method is fast and 
precise enough and can be easily applied for recently proposed enchantments of Balitsky- 
Kovchegov equation where coordinate dependence of gluon forward scattering amplitude 
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interesting by itself, and this enchantments can not be rewritten easily in momentum space. 
First we rewrite Balitsky-Kovchegov equation in terms of variable n(r, x) = 1 — N{r, x) 

It is clearly seen that method used for numerical calculation of integral in can be 
used here too, since integral in (|20j) has very similar structure, however since n(r, x) has 
not singularity at r = we need not exclude area around this point. Moreover since 
n(0, x) = 1 this gives some not vanishing contribution to integral H2U|) we shouldn't do 
that. However we could easily drop region \z\ > z^ax (where Zmax some big value) from 
integration since n(r, x) falls rapidly. Therefore we choose transformation in the following 
form 

(..^i^f^-^i^^ (21, 



,Z — C IT — tpQ ^ 

where 'i/'o = arcos{min(l, (ro/z)) and < ro < r and constants a,b,c will be chosen to 
follow to the following conditions 

— -6 = 1 (22) 
c 

^^6 = 



a 



-b = 0.5 

z' — c 

where z' is some characteristic scale for n{z,y). After this integration can be done by 
adaptive Cubo method. Next we should decide how to do rapidity integration in ()2U() . It 
is oblivious that since n(z, y) will fall with rise of y rises and approach asymptotically to 
zero it is not possible to use explicit integration schema here since it breaks this asymptotic 
behavior easily. Therefore we'll use second order implicit integration method for integration 
since it gives necessary precision in integration and does not violate asymptotic behavior. 



3 Cronin momentum in p + A, d -\- A and A -\- A 
collision 

There are three cases in high energy nuclear collision in which we should check Cronin 
momentum dependency: p + A, d + A, A + A. But first we should define few parameters 
and initial conditions, i.e. form of gluon distribution function for p, d and A and fix value of 



7 




Figure 2: Forward gluon scattering amplitude at different x evolved by Balitsky-Kovchegov 
equation 
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as- For nucleus A (which we suppose to be Au) we'll take unintegrated gluon distribution 
functions in the form (|^. (|1U|) with forward scattering amplitude evolved from initial form 

Noir) = 1 - e-'^'<3.(-)2 in(i/rA+.)/4 (23) 

given at xq = 0.01 with initial saturation momentum Qq^ = 2Gev and regularization 
parameter e = 0.1 (obtained results does not depend on exact value ofeifO<e<l) 
For deuteron we set same initial condition as in (|23|) with saturation momentum rescaled 
according mass number difference. Finally for proton we chose simple 'bremsstrahlung' 
form for unintegrated gluon distribution function. There is only a<j value left to set. Since 
we consider Balitsky-Kovchegov equation with fixed strong coupling we can not set it using 
some scale. However since (jH)),© exhibit same asymptotic behaviour we can relate value 
of A obtained in high energy region and dg with the following relation: 

A = d^ (24) 

7c 

However the only known A = 0.3 was fixed by DIS data at x=0.01 8 . The asymptotic one 
should be lesser that A = 0.3 but scaling law © does not contain all terms and does not 
allow determine exact value of A. Therefore we determine A and consequently a by the 
following method. Fix some some arbitrary value d, calculate forward scattering amplitude 
evolved by Balitsky-Kovchegov equation, determine linear component of scaling law © at 
X = 0.01 and rescale a to fit A = 0.3. After that we can calculate energy and rapidity 
dependency of Cronin momentum. 

Energy dependency of Cronin momentum for p + A collision is shown in figures Ol and 4. 
It is clearly seen that logarithm of Cronin momentum calculated using unintegrated gluon 
distribution function (jUJ higher that Cronin momentum calculated using unintegrated gluon 
distribution function (|10() . The ratio of Cronin momentum for gluon distribution functions 
(jU, ifTU]) is approximately constant and equal exp{0.37). The same also holds also for d + A 
and 74 + ^ collision where Cronin momentum exhibits the same energy behaviour. In A+p 
collision rapidity dependence of Cronin momentum mimics exactly energy dependence since 
gluon distribution function of proton was chosen in simple 'bremsstrahlung' form. 

Rapidity dependence for A + ^ and A + d collisions is shown in figures Q and Q . It is 
clearly seen that Cronin momentum in A + A collision does not have rapidity dependency 
which is quite unexpected behaviour while for A + d we have dependency which consistent 
with (HJ. 
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Figure 3: Energy logarithm dependence of Cronin momentum for A + p collisions at different 
rapidities with cross section calculated using gluon distribution function (jH)) (dotted line) and 
ffTUj) (simple hne) 
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Figure 4: Rapidity dependence of Cronin momentum for A + p collisions at different energies 
with cross section calculated using gluon distribution function ^ (dotted line) and (fTUI) (simple 
line) 
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Figure 5: Rapidity dependence of Cronin momentum for A + A collisions at different energies 
with cross section calculated using gluon distribution function (fTUj) 
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Figure 6: Rapidity dependence of Cronin momentum for A + d collisions at different energies 
with cross section calculated using gluon distribution function (fTUj) 
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4 Conclusion 



Even usage of Balitsky-Kovchegov equation for forward scattering amplitude leads to en- 
ergy behavior of Cronin momentum related with scaling law by dimensional consideration. 
It is clearly seen that although in simple model ^ with scaling law ^ Cronin momen- 
tum calculated using different gluon distribution function have slightly different energy 
behaviour in case of more natural gluon distribution function defined throw solution 
of Balitsky-Kovchegov equation there is no difference between Cronin momentum energy 
behaviour. It is easy to understood since Balitsky-Kovchegov equation noes not contain 
soft scale except in initial condition. 
However: 

• We completely ignore the fact that equation (jSJ can be written using running coupling 
constant. It was shown that this modification changes greatly scaling law © but 
does not break scaling at all. This fact alone should set the same Cronin momentum 
energy behaviour for both definition of unintegrated gluon distribution. Since this 
modification leads to scaling law which does not restores to © even at asymptotically 
high energies it could be easily tested experimentally ether with more precise Cronin 
momentum rapidity dependence or with higher energy measurement. This allows to 
test if a way to introduce running coupling constant is the right one. 

• It was shown in 9 what finite nucleus size sets yet another scale related to its radius 
R. Since this scale is soft it should leads to different results with unintegrated gluon 
distribution functions Q, (|lflj) in both fixed or running coupling case even if we 
consider forward scattering amplitude Nq evolved by Balitsky-Kovchegov equation. 

• In Cronin momentum rapidity behaviour for A + A collisions it is cleanly seen that 
there is no rapidity dependence of Cronin momentum. However since we should 
take into account two different dipole forward scattering amplitude (quark and gluon 
dipole) which obeys different scaling law this behaviour could be easy broken and 
should be tested. 
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